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' We calculate universal finite-size scaling functions for systems with an n-component order pa- 

rameter and algebraically decaying interactions. Just as previously has been found for short-range 
interactions, this leads to a singular £-expansion, where e is the distance to the upper critical di- 
mension. Subsequently, we check the results by numerical simulations of spin models in the same 
, universality class. Our systems offer the essential advantage that e can be varied continuously, al- 

■ lowing an accurate examination of the region where e is small. The numerical calculations turn out 
' to be in striking disagreement with the predicted singularity. 
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In order to analyze numerical results obtained by Monte Carlo or transfer-matrix studies of phase transitions and 
critical phenomena, finite-size scaling ||^ is a very widely used technique. This hypothesis, which was for the first 
time properly formulated by Fisher allows the extrapolation of the properties of finite systems, which do not 
exhibit a phase transition, to the thermodynamic limit. In 1982, Brezin [^| achieved a breakthrough by showing that 
the finite-size scaling laws can actually be derived from renormalization-group (RG) theory, provided that the RG 
, equations are not singular at the fixed point. This implies a breakdown of finite-size scaling for dimensionalities d > 4 
and consequently an expansion of the finite-size scaling functions in powers of e = 4 — d is singular at e = 0. This 
rather surprising result was confirmed by explicit calculations for the n- vector model in the large-n limit. In addition, 
■ it follows from Ref. ^ that the finite-size scaling relation for the free energy is a universal function depending only on 
' two nonuniversal metric factors, without an additional nonuniversal prefactor. This result was derived from different 
I arguments by Privman and Fisher Q and subsequently confirmed analytically for the spherical model |^ . Pioneering 
Q ■ work [||J^ then showed that a field-theoretic calculation of finite-size scaling functions is actually possible. Specifically, 
C^, Brezin and Zinn- Justin |^ developed a systematic e-expansion for these functions. Unlike the standard expansion in 
powers of e for critical exponents and scaling functions of bulk properties, one finds, for a fully finite geometry, an 
7— I \ expansion in powers of ^/e. More recently, Esser et al. |^ introduced a promising perturbation approach at fixed d 
^ ' which is applicable below the critical temperature as well. However, here we focus on the expansion in e and in 
particular on the singular nature of this expansion. 

The systems under consideration have an n-component order parameter cj) with 0(n) symmetry and periodic 
' boundary conditions. A quantity of central interest is the amplitude ratio Q = limi^oo(0L)^/('^i)i which is directly 
I , related to the cumulant introduced by Binder At the critical temperature it takes a universal, although 
• geometry-dependent, value (cf. Ref. In Ref. an expansion for Q{Tc) was obtained in powers of ^/e, up to 

^\ 0{e), which is shown in Fig. |l| for n — 1, along with numerical results for integer d. Given the low order of the 
"Y*^ • expansion and the fact that it can only be checked for integer values of e, hardly any conclusions can be drawn from 
' a comparison to numerical results and any confirmation of the singular nature of the e-expansion will have to wait 
^ until the RG calculation has been carried to substantially higher order. Thus, we propose a different route: namely, 
I ' we replace the short-range (SR) forces by long-range attractive interactions decaying as a power law, J(r) cx 

where < a < 2. It has been shown in Refs. that these systems have an upper critical dimension dc — 2a 

and that for d < dc the critical exponents can be calculated in terms of an expansion in e' = 2a — d, very similar 
^ ' to the original e-expansion, which is recovered for tr — > 2. Since the upper critical dimension now is a continuous 
parameter, we have the opportunity to verify e'-expansion results for arbitrarily small e'. So, even if actual physical 
realizations of this system may be scarce, it constitutes a very valuable mathematical model. Interestingly, the finite- 
size scaling functions for the spherical model with power-law interactions have been calculated by several authors p^ ] 
$H for d/2 < a < d, but the nature of a possible singularity in the limit a — > d/2 appears not to have been examined. The 
first part of this paper is therefore devoted to a generalization of the treatment of Ref. |^ to systems with algebraically 
decaying interactions. For notational convenience we redefine e = e' = 2tT — d. Throughout our analysis we will closely 
adhere to the approach outlined in Ref. 0]. Additional details can also be found in Ref. [p3, Ch. 36. 
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We consider a system with a hypercubic geometry, with hnear dimension L and periodic boundary conditions. It 
is represented by the following Landau-Ginzburg- Wilson Hamiltonian in momentum space, 



1 1 

4! 77 



Uq 



ki,k2,k3 i,j 



(1) 



where the factor fc"^ (0 < cr < 2) arises from the isotropic long-range interactions. Compared to the SR case, it 
leads to the general replacement k"' in all propagators. The indices i,j refer to the components of the field. 

It has been well established that this system belongs to the same universality class as a discrete spin model with 
algebraically decayin g in teractions; in particular the e-expansion results [ p^ , pl| for the critical exponents have been 
confirmed (see Refs. [|l4||l5| references therein). Due to the finite geometry all components of the wave vectors 
are integer multiples of The sums run to infinity, which corresponds to a vanishing lattice spacing a; however, 
the ratio L/^ is finite, whereas L/a and f/a are both sent to infinity Expectation values are computed from a 
partition function in which is replaced by an effective Hamiltonian, consisting of the exactly calculated k = 
(homogeneous) mode contribution and a perturbatively calculated part, which contains the contribution of all nonzero 
modes. Only the latter contribution is affected by the changeover to long-range interactions, cf. Ref. |l^. To one- 
loop order it consists of a shift of the critical temperature and a renormalization of the coupling constant. Higher 
operators do not contribute at this order. We introduce the dimensionless coupling constant go — /i^^uq and work in 
the system of units where /i = 1. The parameter ro is split into roc + t, where t = ro — roc oc T ~ Tc and we require 
t > 0. The RG calculations are carried out in the minimal subtraction scheme [^,0, where it is a crucial ingredient 
of the calculation that, despite the quantization of all momenta, the UV divergences are taken care of by the bulk 
renormalization constants. As we do not go beyond one loop, we have ignored the field renormalization constant. 
The leading contribution to the shift of Tc is 
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where the prime indicates that the k 
be rewritten as 



mode is omitted from the sum. In the Schwinger parametrization this can 
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where — X^iLi ^^"^ '^'^ have omitted the prefactor ^^^i^go- For d> a, the UV divergence in Eq. (^ is reflected 
by the divergence of the integral at small s. Thus we isolate this divergence by rewriting Eq. (^) as 
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which is finite. 5^-1 — 27r''/^/r(|) is the surface area of a d-dimensional unit sphere. The second term in Eq. (^) can 
be continued analytically for 5R(d) > a and has a simple pole for d — 2a (e = 0). Upon expansion around this pole 
we find for Eq. (|) 



2t 



(47r)'^r(cr)£ (47r)^r((T) \<j 



h\t- [In47r + 1'((T)] 



(27rL)' 



:/i(2cT,a,i)+0(e) 



(6) 



where ^'(cr) denotes the Digamma function. The addition of the (/)^-insertion counterterm leads to the replacement of 
t by tZ^2. To one-loop order, the renormalization constant has the same form as for SR interactions |l8|; the pole is 
canceled and the shifted reduced temperature is given by 
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with go = 2[(47r)'^r((T)] ^5o{l + |£[ln47r + ^'(cr)] + ©(e^)}. The renormahzed coupUng constant g is calculated in 
a similar fashion. The leading finite-size contribution is given by — ^^ffo^"'' X^k '(I'^-T + which has a UV 

divergence for di{d) > 2a. The 1/e pole is canceled by the counterterm, where the renormalization constant to 
one-loop order is given by Zg = 1 + ^^^§^5- After some algebra, we find for the rcnormalizcd coupling constant 
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Equations (|^) and ^ suffice to calculate the finite-size scaling functions close to criticality to 0{e). To this order, 
the fixed-point value of 50 only differs from the SR case in the definition oi e, g^ = + 0{e^). We are now able 
to calculate iL'^/'^g~^/'^ at the fixed point, which, as we shall see shortly, is the parameter appearing in the finite-size 
scaling functions. This also provides a simple consistency check for our calculations, since all factors InL have to 
disappear upon introduction of the appropriate scaling variable y = tL^/'^ . We find \/v ~ a — ^7^^ + Oie^), which 
indeed agrees with Ref. , and the final expression is 
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We are particularly interested in the amplitude ratio Q at criticality, for n = 1. The even moments of the magnetization 
distribution are calculated from {{(t>^Y) = d</> (02)p exp[-S'((?!))]/ d0 exp[-5(</.)], with 5(0) = L'^{\t4>^ + 
jfg(j)'^). We carry out the rescaling (f> — > {L'^g)^'^/'^^' and expand in terms of the parameter x = tL'^^'^g^^l'^ . Elementary 
algebra leads to 
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At criticality, y = and x takes the value 
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A comparison to Eq. (3.33) in Ref. ||] shows that xq only differs from the SR case by a redefinition of e, a geometric 
factor, and the integral 1\ and thus we see that the singular structure is preserved in the generalization to long- 
range forces. For completeness we remark that one may also calculate xq by carrying out all manipulations at 
criticality. This permits a different parametrization and leads to the same expression for xq, in which /i(2cr, cr, 0) 
is replaced by /i(2cr, cr)/r(a/2), where h(d,(j) = /o°° duu''/2-i[(X^^_^ e"""')'' - 1 - {-Kluf^ As a side note, 

we conjecture /i(4, 2), which has been evaluated numerically in Refs. mi^,|I^, to be exactly equal to —8 In 2. For 
lower dimensionalities, the upper critical dimension shifts toward smaller values of cr, and numerical evaluation yields: 
/i(d, d/2, 0) = -2.92, -3.900, -4.8227 for d = 1, 2, 3, respectively. Substitution into Eqs. (|l|), (|ll|) suggests that for 
each of these values a reasonable convergence may be expected for e < 1. 

In order to verify these predictions, we have carried out extensive Monte Carlo simulations of spin models with 
n = 1 and algebraically decaying interactions, for d = 1, 2. Accurate results could be obtained by means of an efficient 
cluster algorithm [Q. We investigated system sizes 10 < L < 150 000 for d = 1 and 4 < L < 400 for d = 2, for 
several values of the decay parameter a. These were chosen such that < e ^ 1, where it should be noted that for 
very small e the analysis is hampered by strong corrections to scaling. Simulational details can be found in Ref. [p^ , 
where the classical regime < cr < d/2 is discussed. The numerical results were analyzed using an expression similar 

to Eq. (13) in Ref. IJ, Ql{T) = Q + ritLy* + rat^L^^* H + siLV^ H . Here yt and yi are the thermal and 

leading irrelevant exponent, respectively, and Si are nonuniversal coefficients and the ellipses denote higher-order 
terms. An extensive analysis of the data will be presented elsewhere. The resulting estimates for Q{Tc) are shown 
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in Figs. H, H For the one-dimensional case (Fig. H), the region 0.1 < e < 0.9 has been covered. For a ^ 1, the 
data points approach QiT^) = 1, in agreement with the occurrence of a KosterUtz-Thouless transition at a — 1 [^ . 
However, for < e ^ 0.5 the numerical results are described by a perfectly linear dependence on e, in strong contrast 
with the predicted square-root behavior. This discrepancy is reinforced by the two-dimensional results (Fig. ]^), which 
are also well described by a linear relation for < e < 1.2. For larger cr, the error bars increase, signaling the crossover 
to short-range criticality. In contrast, it should be noted that for the critical exponents good agreement with the RG 
predictions has been reached for at least < e ^ 0.5, both for d — I and d = 2 JTsf . 

In summary, we have calculated universal finite-size scaling functions to second order in ^/e for systems with 
algebraically decaying interactions. These calculations are essentially a generalization of those for systems with short- 
range interactions |^ and exhibit the same singular dependence on e. Subsequently, we have compared our results 
to accurate simulations for one- and two dimensional systems belonging to the same universality class as the field- 
theoretical Hamiltonian. The presence of long-range interactions offers the advantage that e is a continuous parameter, 
so that one can reach the regime where the convergence of the e-expansion has not to be doubted. Nevertheless, no 
agreement is found: the numerical data exhibit a linear rather than a square-root dependence on e. Currently, we do 
not have an explanation for this striking discrepancy. Although higher-order terms might yield some improvement, 
it is difhcult to envisage that this would fully resolve the problems. It would be very remarkable if the apparent 
linear variation over such a wide range in Figs. |^ and ^, which includes the point where our \/e-expansion starts from, 
were accidental. In Ref. p2[ , it has been suggested that Q contains nonuniversal contributions, depending on the 
cutoff used in the integration of the order parameter probability distribution. However, apart from the validity of 
this suggestion, it is difficult to envisage how this would lead to the (dis)appearance of a square-root contribution in 
the e-expansion. Furthermore, it is an open question to what extent the breakdown of the field-theoretic description 
of finite-size scaling for d > 4 ||l^ influences the nature of the e-expansion. We feel that an understanding of these 
problems is of some significance for the understanding of finite-size scaling of critical phenomena in general. 
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Note: After completion of this work, we became aware of Ref. |23| , which contains a different generalization of 
the results of Ref. It allows for the presence of both short-range and long-range interactions, but imposes the 
additional restriction that a = 2 — 0(e). 
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FIG. 1. The amplitude ratio Q{Tc) for systems with short-range interactions. The dashed line shows the -ye-expansion of 
Ref. |. 
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FIG. 2. The critical value of the amplitude ratio Q as a function of the decay parameter a for d = 1. The corresponding 
values of e are shown at the upper axis. The data in the regime < cr < 0.5 are taken from Ref. |14[. 
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FIG. 3. The analog of Fig. | for d = 2. 
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